Studies of Topological Distributions of Inclusive Three— and Four— Jet Events in pp 
Collisions at ^/s = 1800 GeV with the D0 Detector 
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Abstract 

The global topologies of inclusive three- and four-jet events produced in pp interactions are 

described. The three- and four-jet events arc selected from data recorded by the D0 detector at 
the Tcvatron Collider operating at a center-of-mass energy of y's = 1800 GeV. The measured, nor- 
malized distributions of various topological variables are compared with parton-lcvel predictions of 
tree-level QCD calculations. The parton-level QCD calculations are found to be in good agreement 
with the data. The studies also show that the topological distributions of the different subprocesses 
involving different numbers of quarks arc very similar and reproduce the measured distributions well. 
The parton shower Monte Carlo generators provide a less satisfactory description of the topologies 
of the three- and four-jet events. 
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INTRODUCTION 



The Fermilab Tevatron Collider provides a unique opportunity to study the properties of strong interactions in pp 
collisions at short distances. The hard scattering is described by the theory of perturbative Quantum Chromodynamics 
(QCD) |^-|^ and has been studied extensively in the last decade |^,p|. Within the context of QCD, the hard process 
is described as a point-like scattering between constituent partons (quarks and gluons) of protons and anti-protons. 
The scattering cross sections can be written in expansions in terms of powers of the strong coupling constant as 
convoluted with parton momentum distributions inside the nucleon. The lowest order term corresponds to the 
production of two-parton final states. Terms of order and in the expansion imply the existence of three- 
and four-parton final states, respectively. Colored partons from the hard scattering evolve via soft quark and gluon 
radiation and hadronization processes to form observable colorless hadrons, which appear in the detector as localized 
energy deposits identified as jets. High energy jets originating from partons in the initial hard scattering process are 
typically isolated from other collision products. They are expected to preserve the energy and direction of the initial 
partons, and therefore the topologies of the final jet system are assumed to be directly related to those of the initial 
parton system. 

The cross section and angular distributions for two-jet events have been successfully compared with the predictions 
of QCD 1^,1). A study of three- and four-jet events allows a test of the validity of the QCD calculations to higher 
order (a^ or beyond) and a probe of the underlying QCD dynamics. This paper explores the topological distributions 
of three- and four-jet events. The distributions provide sensitive tests of the QCD matrix element calculations. 
Topological distributions for the three- and four-jet events have been published previously by the UAl, UA2 and CDF 
Collaborations However, all of these studies imposed requirements on the topological variables themselves, and 

therefore significantly reduced the phase space under study. This paper extends these studies to previously untested 
regions of phase space for a large number of topological variables. The measured normalized distributions, without 
restrictions on the topological variables themselves, are compared with the QCD tree-level matrix element calculations. 
The predictions from simple phase-space matrix elements are shown as a comparison, and the distributions of QCD 
subprocesses involving different numbers of quarks are also examined. Finally, the data are compared with the 
predictions of three parton shower event generators. 

DEFINITION OF TOPOLOGICAL VARIABLES 

The topological variables used in this paper are defined in the parton or jet center-of-mass system (CMS). The 
definitions refer to partons and jets interchangeably. The partons are assumed to be massless and the jet masses are 
ignored by using the measured jet energies as the magnitudes of jet momenta. 

The topological properties of the three-parton final state in the center-of-mass system can be described in terms of 
six variables. Three of the variables reflect partition of the CMS energy among the three final-state partons. The other 
three variables define the spatial orientation of the planes containing the three partons. It is convenient to introduce 
the notation 1 + 2^3 + 4 + 5 for the three-parton process. Here, numbers 1 and 2 refer to incoming partons while 
the numbers 3, 4 and 5 label the outgoing partons, ordered in descending CMS energies, i.e. E^> E^> E^. The final 
state parton energy is an obvious choice for the topological variables for the three-parton final state. For simplicity, 
Ei {i = 3,4, 5) is often replaced by the scaled variable Xi {i = 3, 4, 5), which is defined by Xi = 2Ei/^/§, where Vi is 
the center-of-mass energy of the hard scattering process. By definition, 2:3 + 2:4 + 2:5 = 2. The scaled parton energies 
Xi and the angles between partons {ujjk, j, A: = 3, 4, 5) for the three-parton final state have the following relationship: 

(1) 



sm UJ34 + sm LUi5 + sm 



where i,j,k = 3,4,5 and i ^ j ^ k. Clearly, the internal structure of the three-parton final state is completely 
determined by any two scaled parton energies. The angles that fix the event orientation can be chosen to be: (1) the 
cosine]^ of the polar angle with respect to the beam (cos 6*3) of parton 3, (2) the azimuthal angle of parton 3 (^3), and 
(3) the angle between the plane containing partons 1 and 3 and the plane containing partons 4 and 5 ("0*) defined by: 

cosV.* = %^%^i^, (2) 



^Unless otherwise specified, the absolute values of the cosines of polar angles are implied throughout this paper. 
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where pi is the parton momentum. Figure |i| illustrates the definition of the topological variables for the three- 
parton final state. For unpolarized beams (as at the Tevatron), the ^3 distribution is uniform. Therefore, only four 
independent kinematic variables are needed to describe the topological properties of the three-parton final state. In 
this paper, they are chosen to be X3, x^, cos^g and "0*. 

Another set of interesting variables is the scaled invariant mass of jet pairs: 



^ = J XiXj {1- COS uJij)/ 2 i,j = 3,4,5 and i 7^ (3) 
V s * 

where my is the invariant mass of partons i and j and Uij is the opening angle between the two partons. The 
scaled invariant mass (/Zy ) is sensitive to the scaled energies of the two partons, the angle between the two partons, 
and the correlations between these variables. Using dimensionless variables and making comparisons of normalized 
distributions minimizes the systematic errors due to detector resolution and jet energy scale uncertainty and therefore 
allows a direct comparison between data and theoretical calculation. 

The four-parton final state is more complicated. Apart from the CMS energy, eight independent parameters are 
needed to completely define a four-parton final state in its center-of-mass system. Two of these define the overall event 
orientation while the other six fix the internal structure of the four-parton system. In contrast to the three-parton 
final state, there is no simple relationship between the scaled parton energies and the opening angles between partons. 
Consequently, the choice of topological variables is less obvious in this case. In this paper, variables are defined in a 
way similar to those investigated for the three-parton final state. The four partons are ordered in descending CMS 
energy and labeled from 3 to 6. The variables include the scaled energies (xi, z = 3, 6), the cosines of polar angles 
(cos^*, i — 3,..., 6) of the four jets, the cosines of the opening angles between partons (coswy, i,j = 3,.. .,6, and 
i 7^ j), and the scaled masses (^^ = rriij/y/i, i,j = 3, ...,6 and i 7^ j) of parton pairs. In addition, two variables 
characterizing the orientation of event planes are investigated. One of the two variables is the "Bengtsson-Zerwas" 
angle (xbz) pi[ defined as the angle between the plane containing the two leading jets and the plane containing the 
two non-leading jets: 

{P3 X Pa) ■ {P5 X pe) 

coaxBZ^-rz ^77^ (4) 

\P3 X pa\\p5 X Pel 

The other variable is the cosine of the "Nachtmann-Reiter" angle {cos 9nr) (is) defined as the angle between the 
momentum vector differences of the two leading jets and the two non-leading jets: 

cos 9 NR^ = — i~: ^TTT^ ^-T- (5) 

\P3 ~Pi\\P5~P6\ 

Figure H illustrates the definitions of xbz and Onr variables. Historically, xbz and cosOnr were proposed for e+e" 
collisions to study gluon self-coupling. Their interpretation in pp collisions is more complicated, but the variables can 
be used as a tool for studying the internal structure of the four-jet events. 



THE THEORETICAL MODEL 



The cross section for the production of the n-parton final state 1 + 2 ^ 3 + • ■ ■ + (n + 2), in pp collisions at a 
center-of-mass energy ^/s is described by the following expression: 

a„ = ^ /" fi{xi)f^{x2)\M^\''<i>ndxidx2, (6) 

where the sum runs over all possible 1 + 2 ^ n— parton subprocesses. The functions fi{xi) and f2{x2) are the parton 
density functions of the incoming partons, |M"p represents the matrix elements of the subprocess, and <&„ is the 
n-body phase space. Theoretically |M"p is well behaved if calculated to all orders in the expansion. At present, 
this calculation is technically not possible and one has to deal with truncated expansion. As a result, diverges 
when the energy of any final state parton or the angle between any two partons approaches zero. The singularities 
in |M"p cause poles in the topological distributions. In comparison, a phase-space model in which |M"p cx l/s"~^ 
where s = X1X2S does not have singularities in the matrix element, therefore, the topology of the model is determined 
by the phase space $«. In this paper, the distributions from the phase-space model are used as references for the 
comparisons between the data and QCD. 
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Presently two approaches for modeling perturbative QCD for multi-jet production exist. The straight-forward 
method is the matrix element method, in which Feynman diagrams are calculated order-by-order in as- Technical 
difhculties have limited the calculations to the tree-level of the relevant processes. The exact tree-level matrix 
element calculation for the three-parton final state has been available for some time (l3|. The complete tree-level 
matrix element calculations for up to five final state partons have been recently calculated by Berends, Giele and 
Kuijf (BGK) [Q using a Monte Carlo method. The other commonly used approximate calculations are those of 
Kunszt and Stirling (KS) and of Maxwell The perturbative QCD calculations have been incorporated into 
several partonic event generators. The exact tree-level matrix elements calculations for up to five jets are implemented 
in the NJETS program. PAPAGENO jl7| implements an exact matrix element calculation of tree-level contributions 
for final states with up to three partons and provides KS and Maxwell approximations for up to six partons. These 
approximations are used in part to speed up the calculations, in view of the complicated exact matrix elements. For the 
analysis described in this paper, the NJETSprogram is used to calculate QCD predictions while the PAPAGENOprogram 
is used as a cross check and to calculate distributions from the phase-space model. 

The second approach is based on the parton shower scheme. In this method, the hard scattering begins with two 
initial outgoing partons. An arbitrary number of partons are then branched off from the two outgoing partons and 
the two incoming partons (backward evolution) to yield a description for multi-parton production, with no explicit 
upper limit on the number of partons involved. The parton shower picture is derived within the framework of the 
Leading Logarithmic Approximation (LLA) ]T^ . As a result of the approximation, the reliability of the parton 
shower is expected to decrease as parton multiplicity increases. Many parton shower Monte Carlo event generators 
are available. In this paper, HERWIG 5.8 ISAJET 7.13 and PYTHIA 5.6 are used. 



THE DATA SAMPLE 



The data used in this analysis were collected with the D0 detector during the 1992-1993 Tevatron run at a 
center-of-mass energy of 1800 GeV. The D0 detector consists of a central tracking system, a calorimeter, and muon 
chambers. Jets are measured in the calorimeter, which has a transverse segmentation of Arj x = 0.1 x 0.1. The 
jet energy resolution is typically 15% at Et=50 GeV and 7% at £'t=150 GeV (HJ. The jet direction is measured 
with a resolution of 0.05 in both rj and cj). With the hermetic and uniform rapidity coverage (—4.5 < rj < 4.5) of the 
calorimeter, the D0 detector is well suited for studying multi-jet physics. A detailed description of the D0 detector 
can be found elsewhere |^ . 

The events used in this study passed hardware (Level 1) and software (Level 2) energy-cluster based triggers. In 
addition, a Level hardware trigger required that vertices along the beam line be within 10.5 cm of z = 0. The Level 1 
trigger was based on energy deposited in calorimeter towers of size A77 x Acj) = 0.2 x 0.2. The events were required to 
have at least two such towers with transverse energy (Et) above 7 GeV. The successful candidates were passed to the 
Level 2 trigger, which summed transverse energies of calorimeter towers in a cone of radius 7e(= (Ar;2+A(/,2)i/2) = o.7. 
The Level 2 trigger selected those events with at least one such cone, built around the Level 1 trigger tower, with 
transverse energy above 50 GeV. The total effective luminosity used in this analysis is 1.2 pb~^. The trigger efficiency 
for events with at least one jet with Et > 60 GeV is above 90% ^l|. A detailed description of the trigger can be 
found elsewhere p^ . 

The offline reconstruction uses a fixed-cone jet algorithm with TZ — 0.7, similar to the algorithm used in the Level 2 
trigger. The jet reconstruction begins with seed calorimeter towers of size A77 x A(f) — 0.1 x 0.1 containing more than 
1 GeV transverse energy. Towers are represented by massless four-momentum vectors with directions given by the 
tower positions and event vertices. The four momenta of towers in the cone around the seed tower are summed to 
form the four-momentum vector of the jet. The jet direction is then recalculated using tower directions weighted by 
their transverse energies. The procedure is repeated until the jet axis converges. For two overlapping jets, if either 
jet shares more than 50% of its transverse energy with the other jet, the two jets are merged. Otherwise they are 
split and the shared transverse energy is equally divided between the two jets. The final jet Et is the sum of the 
transverse energies of towers within the cone, while the jet direction is determined by the jet four-momentum vector 

{E,E^,Ey,E^), i.e., e = cos~^{E,J^El + E^ + El), = iim-^{Ey/E^) and -q = - lntan(6'/2). 

The jet energy scale has been calibrated using direct photon candidates by balancing jet Et against that of the 
photon candidate. The electromagnetic energy scale was determined by comparing the measured electron pair mass 
of Z — > e^e~ events with the Z mass [25| measured by e^e~ experiments. The calibration takes into account the 
effects of out-of-cone particle showering using shower profiles from test beam data as well as the underlying event 
using events from minimum-bias triggers. Details can be found in Ref. p^ . 

After energy corrections, jets are required to have Et greater than 20 GeV and lie within a pseudorapidity range 
of —3.0 to 3.0. The pseudorapidity is calculated with respect to the event vertex determined from tracks measured 
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by the central tracking detector. Jets passing the above criteria are ordered in decreasing Et- The Et of the leading 
jet must be greater than 60 GeV to reduce possible trigger bias and threshold effects. 

Three-jet events are selected by further demanding that there be at least three jets. This leaves about 94,000 
events in the sample. The separation ATZ between jets is required to be greater than 1.4, which is twice the cone 
size used, to avoid systematic uncertainty associated with the merging/splitting of the cone jet algorithm. This 
requirement removes events with overlapping jets and therefore ensures good jet energy and direction measurements. 
Approximately 70% of the events pass this requirement. The invariant mass distribution of the three highest Et jets 
is shown in Fig ^. Also shown is the distribution from the exact tree-level calculations of perturbative QCD. The 
overall agreement between the data and QCD distributions is good with the exception of the low mass region, where 
the threshold and resolution effects are important. To reduce possible bias in this region, the invariant mass of the 
three leading jets is required to be above 200 GcV/c^. After all selection criteria, a sample of about 46 thousands 
three-jet events remains. The surviving events are then transformed to the CMS frame of the three leading jets. Any 
other jets in the event are ignored. The jets are re-ordered in descending energy in their CMS system. The topological 
variables (xa, x^, cos 9^ and ^p*) are calculated. Unlike previous studies by other experiments, no requirements on 
these topological variables are imposed. If the topological requirements similar to those in Ref. were imposed, the 
three-jet event sample would be reduced by more than a factor of ten. 

Four-jet events are selected in a similar manner. Events are required to have at least four jets, which results in 
a data sample of 19,000 events. The ATZ between any jet pair is required to be greater than 1.4, reducing the data 
sample to about 8,400 events. As in the selection of the three-jet events, the invariant mass of the leading four jets 
must be above 200 GeV/c^. The mass distribution before this requirement is applied is shown in Fig. |[ A total 
of 8,100 events remains in the four-jet event sample. The four leading jets of the remaining events are boosted to 
their center-of-mass system, and are ordered in decreasing energy. Additional jets, if present, are ignored. The 
topological variables are calculated using the four boosted momentum vectors after ordering in decreasing energy. No 
requirements on the topological variables are imposed. 



PREDICTIONS OF THEORETICAL MODELS 



The partonic event generator NJETS is used to calculate the exact tree-level QCD distributions. The PAPAGENO 
program is used to calculate the distributions of the phase-space model as well as the approximate calculations of KS. 
Unless otherwise specified, the parton distribution function used in the calculations is MRS (bcdms fit) for both 
NJETS and PAPAGENO. The QCD scale parameter is set to 200 MeV and the renormalization scales are set to the 
average transverse momentum of the outgoing partons for both matrix elements and parton distribution functions. 
The outgoing partons are analyzed as if they were observed jets and the selection criteria described above are applied 
to select three- and four-jet events. 

To study the sensitivity to the choice of parton distribution function, the topological distributions of QCD cal- 
culations with different parton distribution functions are compared. For NJETS, the comparisons are made between 
MRS and EHLqJ2^ parton distribution functions. For PAPAGENO the parton distribution functions of MRS |2^] 
and Morfin-Tung ^| are employed. Although the total three- and four-jet cross sections vary by as much as 30% for 
different parton distribution functions, the normalized topological distributions are found to be very insensitive to the 
choice. A typical difference of less than 3% is found for the variables studied. The dependences on the renormaliza- 
tion scale are investigated using the PAPAGENO program. The distributions for the renormalization scales of (1) the 
average transverse momentum, (2) one half the average value of transverse momentum and (3) the total transverse 
energy are compared. Despite large differences (as much as 60%) in the total production cross sections, the differences 
between normalized distributions are very small, typically less than 3%. Combining the effects described above, the 
uncertainty on the theoretical predictions is 4%. 

The fragmentation effect is investigated using the HERWIG event generator. The parton-level distributions for 
three- and four-jet events are compared with the distributions at particle level. For parton shower Monte Carlo 
programs, partons are defined as those quarks and gluons after the parton showering and before the fragmentation. 
The differences between the distributions before and after fragmentation are found to be small, typically at less 
than 4% level. The small fragmentation effect combined with a small detector effect discussed below enable direct 
comparisons between data and theoretical parton-level calculations. 

Both NJETS and PAPAGENO incorporate tree-level calculations for three- and four-parton final states. The effect 
on the normalized distributions due to higher-order loop corrections is expected [po| to be small in the phase-space 
region relevant to the analyses described in this paper. Although both Monte Carlo programs generate exclusive 
events, the three or four jets of the generated events predict the behavior of the leading three or four jets of an 
inclusive data sample |30|. Therefore the data distributions based on the inclusive samples are compared with QCD 
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calculations from exclusive final states in this paper. 



UNCERTAINTIES OF THE MEASURED TOPOLOGICAL DISTRIBUTIONS 

The measured distributions of topological variables are affected by: (a) the trigger efficiency, (b) the detector 
acceptance and resolution, and (c) the uncertainty of the energy scale. However, most of these corrections and their 
uncertainties are minimized by normalizing the distributions to unit area and by selection requirements. In the 
following, residual uncertainties are discussed. 

The non-uniformity of the detector acceptance and of the trigger efficiency in the topological variables and the 
detector energy resolution and angular resolution have direct effects on the measured distributions. These effects are 
estimated using a fast detector simulation program |3l[] which takes into account the detector energy and angular 
resolution and the trigger efficiency as functions of the pscudorapidity and the transverse energy of jets. The bin-by- 
bin correction factors are flat within 5%. 

By definition, the topological variables have a weak dependence on the energy scale since only the scaled energies 
and directions of the jets are used. However, the event selection criteria, such as Et and invariant mass requirements, 
are vulnerable to the energy scale error. The possible distortion of the measured topological variables due to the 
uncertainty in the energy scale is studied by varying the energy calibration constants within their nominal errors. 
The selection procedure described above is repeated for the events calibrated with these modified constants. Apart 
from some low statistics bins, the variations in the measured topological variables are very small. We conservatively 
assign a 3% systematic error on the topological distributions due to energy scale uncertainty. The small variation is 
in part due to the fact that the topological distributions change slowly with the jet Et and the invariant mass of the 
jet system. 

In principle, the measured distributions have to be corrected for detector effects before the data can be compared 
with the theoretical calculations. However, adding the above systematic effects in quadrature, we get a 6% uncertainty 
on the measured distributions. The small detector effects suggest that the data distributions can be directly compared 
with the parton-level distributions of perturbative QCD calculations. In the following, the measured distributions 
with a 6% estimated total systematic error are directly compared with the QCD tree-level calculations at the parton 
level. Finally, we note that changing the jet separation requirement ATI from 1.4 to 1.0 docs not change the degree 
of agreement between the data and QCD calculations. 

THE TOPOLOGIES OF THREE-JET EVENTS 

Figure ^ and Table | show the measured X3 and X5 distributions for the final selection of three-jet events. The three 
jets are labeled in order of decreasing energy in their CMS frame. The average values of and are 0.88 and 0.39 
respectively. The data are compared with the predicted distributions of the exact QCD tree-level calculations and 
the expectations from the phase-space model. The QCD calculations reproduce the measured distributions well for 
the entire range. Unlike the predictions of the phase-space model, the data heavily populate the high 2:3 region and 
have significant contributions at low X5 values, a characteristic of gluon radiation. The decrease in X3 distributions at 
high X3 values is due to the ATZ requirement in the event selection. The bottom plot shows the fractional difference 
between the data and the QCD predictions with dotted lines indicating the estimated 6% systematic error on the 
measurement. The differences between the data and the predictions are generally within the systematic error band. 
The root mean square (RMS) of the fractional differences between the data and the QCD predictions are 3.4% for X3 
and 3.9% for X5. 

The cos ^3 distribution is shown in Fig. ^(a). As in the angular distribution of two-jet events, an angular dependence 
characteristic of Rutherford i-channel scattering is unmistakable. The large angular coverage of the D0 calorimeter 
allows the analysis to cover the entire cos ^3 range, extending the study into a previously untested region of phase space. 
As is evident in the figure, the data are well reproduced by the predictions of the exact QCD tree-level calculations 
over the entire range of cosflg, with a RMS 4.0% of the fractional differences. The phase-space distribution is mostly 
fiat with high cos 9^ bins suppressed as a result of the pscudorapidity requirement in the event selection. The depletion 
in the data and the QCD calculations is compensated by a large cross section in this region and therefore is less visible. 
The measured distribution is shown in Fig. ^(b) together with the results of the exact QCD tree-level calculation 
and of the phase-space model. The phase-space distribution shows depletions at small and large ip* angles, an effect 
of the event selection. However, the data and the QCD distributions are enhanced in these regions because of initial- 
state radiation in which one of the two non-leading jets is close to the beam line. As in the case of the 3:3, X5 and 
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COS ^3 distributions, the overall agreement between data and the QCD tree-level calculations is very good. The RMS 
of the fractional differences is 4.2%. 

The scaled mass distributions are sensitive to the jet energies, the opening angles between jets, and the correlations 
between these quantities. The measured /^34, fi^^ and /i45 distributions for the three-jet event sample are compared 
with the exact QCD predictions in Fig. |^. The QCD predictions agree with the data well, while the differences between 
the data and the phase-space model are large. We also note that some systematic shift in fi^^ and /i45 distributions 
is clearly visible. The RMS's of the fractional differences between the data and the QCD calculations are 3.6%, 6.7% 
and 6.9% for /X35 and 1145 respectively. 

Finally, we note that the KS approximate QCD calculations are essentially identical to the exact tree-level QCD 
calculations for the topological variables studied above. This implies that the topological distributions are insensitive 
to the approximation made in the KS calculations. 



THE TOPOLOGIES OF FOUR-JET EVENTS 



The four measured energy fractions of four-jet events are shown in Fig. |7| and also listed in Table IV. The four jets 
are labeled in order of decreasing energy in their center-of-mass system. Although four scaled energy variables are 
shown, only three of these are independent. The other is fixed by the condition Xi = 2. The measured mean values 
of the four energy fractions are 0.76, 0.61, 0.39, and 0.24. The QCD predictions of the exact tree-level calculations 
are represented by the solid curves and are in an excellent agreement with the data for all four variables. As in the 
three-jet case, the distributions from the phase-space model do not reproduce the data. The fractional differences 
between the data and QCD arc very similar to those of the three-jet events and are not shown for simplicity. 

The cosines of the four polar angles of the four-jet events in their center-of-mass system are compared with 
QCD calculations in Fig. || for the entire range. While the two leading jets tend to be in the forward region, the 
cosine distribution of the least energetic jet is essentially flat, because the jet separation requirement in the event 
selection favors events with other jets in the central region. Although small differences between the data and the 
QCD calculations are visible, the overall agreement is good. Despite the large differences between the data and the 
phase-space model in cos 6*3 and cos 91 distributions, the differences in the other two distributions are relatively small. 

The internal event structure can be further understood by examining the opening angles between jets. Figure ^ 
shows the distributions of the space angle between all possible jet pairs of the four-jet events in their center-of-mass 
system. While the two leading jets are mostly back-to-back, the angles between other jet pairs are distributed widely. 
The depletion in the regions where coswy — > 1.0 is again due to the ATZ requirement in the event selection. The 
structures of the data distributions are well described by the QCD predictions. 

Figure |l^ shows the scaled mass distributions of jet pairs of the four-jet events for both data and the QCD 
calculations. The average scaled mass is 0.65 for the two leading jets and is 0.23 for the two non-leading jets. The 
QCD calculations agree with the data well. Distributions of the phase-space model are generally too narrow and fail 
to reproduce the data distributions. 

Figure pi] compares the measured xbz and cos 9nr distributions with the predictions of the exact tree-level QCD 
calculations as well as those from the phase-space model. The agreement between the data and QCD is generally good 
and the differences between the data and the phase-space model are large. Although the jet separation requirement 
in the event selection favors large xbz-, the data and the QCD distributions have significant contributions in the 
small Xbz region, which corresponds to a planar topology of the four jets. In contrast, the phase-space distribution 
is highly suppressed in this region. The cos On b. distributions for the data and QCD are essentially flat while the 
phase-space model peaks strongly as cos 9 nb. approaches zero. 

For the four-jet events as was the case for the three-jet events, the normalized distributions from the KS approximate 
QCD calculations agree well with the data. 



COMPARISON OF QCD SUBPROCESSES 



At the parton level, five and six partons (including the two initial partons) are involved in the three- and four-jet 
processes respectively. It is difficult, if not impossible, to label quark or gluon jets in the data. However, with Monte 
Carlo event generators, the three-jet cross section can be broken into three subprocesses involving different numbers of 
quarks among the initial- or final-state partons: (1) 0-quark, (2) 2-quark and (3) 4-quark. The predicted fractional 
contributions by NJETS to the total three-jet cross section for the selection criteria described above are 32.9%, 50.8% 
and 16.2% for 0-quark, 2-quark and 4-quark subprocesses respectively. Similarly, the four-jet process can be divided 
into (1) 0-quark (29.4%), (2) 2-quark (49.6%), (3) 4-quark (20.2%) and (4) 6-quark (0.7%) subprocesses. 
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The studies described above show that the QCD calculations agree well with the data. It is therefore interesting to 
examine the topological distributions of these subprocesses. Figures |l^ (a) and (b) show the X4 and cos 91 distributions 
of the three-jet events and Figs. |l2| (c) and (d) show the X5 and cos^g distributions of the four-jet events predicted 
by the exact tree-level QCD calculations (full QCD) and by the QCD calculations of the three subprocesses. The 
full QCD is normalized to unit area and the subprocesses arc normalized to the fractional contribution to their 
respective total cross section. The data distributions are normalized to the respective QCD distributions. The 
distributions of the subprocesses are remarkably similar and agree well with the data. The 6-quark subprocess 
contributes less than 1% of the total cross section of the four-jet events and therefore is not shown in Figs, [ij (c) 
and (d). Nevertheless, the normalized distributions are very similar to those of the other three subprocesses. The 
similarity of the subprocesses is observed in all other variables of the three- and four-jet events investigated in this 
paper. This suggests that the distributions are insensitive to the relative contributions of these subprocesses to the 
total cross section and therefore have weak dependences on the quark-giuon content in parton distribution functions. 
Futhermore, Rutherford characteristics are visible in cos 9* distributions for all subprocesses, implying that the matrix 
elements of these subprocesses are dominated by the t-channel exchange. 

COMPARISON WITH PARTON SHOWER EVENT GENERATORS 

As discussed above, the measured topological distributions of three- and four-jet events are reproduced well by 
the exact tree-level QCD calculations. However, in many investigations, parton shower Monte Carlo event generators 
are used to model multi-jet production. Therefore, it is interesting to compare the data distributions with those 
predicted by parton shower event generators. 

As an example, the X3 and cos ^3 distributions of three-jet events and ^134 and fi^e distributions of four-jet events 
are shown in Fig. ^for the data and for the HERWIG 5.8, ISAJET 7.13 and PYTHIA 5.6 parton shower event generators^. 
The Monte Carlo distributions are calculated using parton jets which are formed by quarks and gluons after parton 
showering and before hadronization. The parton jets are initially reconstructed using a cone jet algorithm implemented 
in the PYTHIA 5.6 program Then the jet direction is redefined using a D0 jet direction definition discussed above. 
Although the parton shower generators describe the general structures of these variables well, differences in details 
are clearly visible. The largest difference is seen in the cos 6*3 distribution. All three parton shower event generators 
show excessive contributions in the forward region. 

To generate three- and four-jet events using the parton shower generators, one has to start with 2^2 processes 
with a pt cut and select events with hard gluon radiation. We note that a large fraction of the Monte Carlo events 
in the forward region which pass the 60 GeV leading jet Et requirement have 2^2 process with pT < 50 GeV/c. 
Presumably the leading jets of the these events are from hard initial-state radiation. This observation suggests that 
the initial-state radiation is not well modeled by these parton shower generators in the phase-space region studied in 
this paper. 

Although only four topological distributions are shown here, we have compared all other variables investigated in 
this paper. Apart from the cos 9* distributions, the HERWIG event generator provides a reasonably good description 
of the data while the differences between the data and the predictions of ISAJET and PYTHIA event generators are 
large in many distributions. Overall, the HERWIG event generator describes the data better than the ISAJET and the 
PYTHIA do. 



SUMMARY 

From the data sample recorded by the D0 detector in pp collisions at y/s — 1800 GeV at the Tevatron during the 
1992-1993 running period, high statistics three-jet and four-jet event samples have been selected. A large number 
of distributions characterizing the global structures of the inclusive three- and four-jet events have been compared 
with QCD calculations of the exact tree-level matrix elements and with calculations of QCD subprocesses involving 
different numbers of quarks. This paper extends earlier studies to previously untested regions of phase space for 
a large number of topological variables. For example, compared with an earlier study |^ of the three-jet events, 
the cos6'3 region studied has been expanded from 0.8 to 1.0, the 2:3 uplimit from 0.9 to 1.0 and the ^/;* range from 



^All parton shower events are generated with a pT = 10 GeV/c cutoff for the initial 2 — > 2 hard scattering, using their default 
parameters. 



9 



20° < tp* < 160° to 0° < tp* < 180° for a minimum thrcc-jct invariant mass of 200 GcV/c^. All comparisons have 
been made with the parton-lcvcl distributions and based on normalized distributions rather than cross sections. 

For the three-jet events, the investigated topological variables arc: the energy fractions carried by the two leading 
jets, the cosine of the leading jet polar angle, the angle between the plane containing the leading jet and the beam 
line, the plane containing the two non-leading jets, and the scaled invariant masses of the jet pairs. In the case of the 
four-jet events, the energy fractions and the cosines of the polar angles of all four jets, the six opening angles, scaled 
invariant masses of jet pairs, and the angles between jet planes have been studied. 

Studies show that the measured topological distributions of the three- and four-jet events are well reproduced by 
the exact tree level matrix elements QCD calculations. The good agreement implies that the topological distributions 
of the three- and four-jet events are determined by the tree-level diagrams and therefore the topological distributions 
are not very sensitive to higher-order corrections. Futhermore, the distributions are found to be insensitive to the 
uncertainties in parton distribution functions and to the quark/gluon flavor of the underlying partons. The dominance 
of the i-channel gluon exchange to a large extent determines the structure of the event. The differences between the 
data and the phase space model are large for most of the distributions. The successful direct comparison between the 
data and the QCD calculations at the parton level reaffirms the assumption that jets closely follow their underlying 
partons at high energies. Finally, we note that apart from the cos 0* distributions, the HERWIG 5.8 event generator 
provides a good description of the measured distributions while the differences between the data and the predictions 
of the ISAJET 7.13 and the PYTHIA 5.6 event generators are relatively large in many distributions. 
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FIG. 1. Illustration of the topological variables: Xi, uiij, 63 and ijj* for the three-jet events. 




FIG. 2. Illustration of xbz and Onr definitions for the four-jet events. 
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FIG. 3. The mass distributions of the selected three- and four-jet events before the mass requirement. The QCD matrix 
element calculations are normalized to the data. 
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FIG. 4. The scaled energy distributions of (a) xz and (b) for three-jet events in their center-of-mass system. Only 
statistical errors are shown. The bottom plot shows the fractional difference of the data from the exact tree-level QCD 
calculation. The RMS of the fractional differences is 3.4% for x^ and 3.9% for xz,. The dotted lines show the estimated 6% 
systematic uncertainty on the measurement. The data are also listed in Table H 
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FIG. 5. The distributions for the three-jet events of (a) the cosine of the leading jet polar angle and (b) the angle t/i* (defined 
in the text) in their center-of-mass system. The RMS of the fractional differences is 3.4% for cos 6'^ and 4.2% for ■(/)*. The 
dotted lines show the estimated 6% systematic uncertainty on the measurement. The data are also listed in Table 
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FIG. 6. The scaled mass distributions of (a) /i34, (b) /X35 and (c) ^45 for the three-jet events in their center-of-mass system. 
The bottom plots show fractional differences between the data and QCD. The RMS's of the fractional differences are 3.6%, 
6.7% and 6.9% respectively. The dotted lines show the 6% systematic uncertainty on the measurement. The data are also 
hsted in Table 
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FIG. 7. The jet energy fraction distributions for the four-jet events in their center-of-mass system. Only statistical errors 



are shown. The RMS's of the fractional differences between the data and QCD are 6.5%, 5.5%, 6.2% and 6.5% for x^, 



and x% respectively. The estimated systematic uncertainty on the measurement is 6%. The data are also listed in Table IV 
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FIG. 8. The distributions of jet polar angle for the four-jet events in their center-of-mass system. Only statistical errors are 
shown. The RMS's of the fractional differences between the data and QCD are 6.8%, 6.2%, 7.5% and 7.9% for cos^l, cos^l, 
cos and cos 65 respectively. The estimated systematic uncertainty on the measurement is 6%. The data are also listed in 



Table 
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cosco 

FIG. 9. Distributions of the space angle between jet pairs for the four-jet events in their center-of-mass system. Only 
statistical errors are shown. The RMS's of the fractional differences between the data and QCD for the six space angles are 
5.4%, 5.6%, 6.4%, 4.7%, 7.6% and 8.0%. The estimated systematic uncertainty on the measurement is 6%. The data are also 
hsted in Table 
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FIG. 10. Distributions of scaled jet pair mass for the four-jet events in their center-of-mass system. Only statistical errors 
are shown. The RMS's of the fractional differences between the data and QCD for the six masses are 7.8%, 7.1%, 6. 8%, 6.3%, 
4.1% and 3.7%. The estimated systematic uncertainty on the measurement is 6%. The data are also listed in Table VII. 
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FIG. 11. The distribution of the angle between the planes of (a) the two leading jets and the two non-leading jets and (b) 
the momentum vector differences of the two leading jets and the two non-leading jets, for four-jet events in their center-of-mass 
system. Only statistical errors are shown. The RMS of the fractional differences is 5.7% for x bz a nd 4.1% for cos 9nr- The 
estimated systematic uncertainty on the measurement is 6%. The data are also listed in Table VIIl. 
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FIG. 12. The (a) X4 and (b) cosOX distributions for three-jet events and the (c) xs and (d) cosflj distributions for four-jet 

events in their center of mass system. The QCD subprocesses are normalized to their fractional contributions to the respective 
total cross section for the selection criteria described in the text. The data are scaled to the normalization of the respective 
subprocess of the QCD calculation. Therefore, only the shapes of the subprocesses are compared. The estimated uncertainty 
on the data is less than 6%. 
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FIG. 13. Comparisons between the data, exact trcc-lcvcl QCD calculations, and HERWIG, ISAJET and PYTHIA Monte 
Carlo predictions. Shown are (a) the scaled energy of the leading jet and (b) the cosine of the leading jet for thrce-jet events, 
and the scaled invariant mass distributions of (c) the two leading jets and (d) the two non-leading jets for four-jet events. 
Only statistical errors are shown. The estimated systematic uncertainty on the measured distributions is less than 6%. 



21 



X3 


l/NdN/dx3 




l/NdN/dxi 


0.66-0.68 


0.069±0.009 


0.075-0.100 


0.155±0.012 


0.68-0.70 


0.409±0.021 


0.100-0.125 


0.417±0.019 


0.70-0.72 


0.751±0.028 


0.125-0.150 


0.680±0.024 


0.72-0.74 


0.997±0.033 


0.150-0.175 


1.125±0.031 


0.74-0.76 


1.480±0.040 


0.175-0.200 


1.488±0.036 


0.76-0.78 


1.828±0.044 


0.200-0.225 


1.881±0.040 


0.78-0.80 


2.190±0.049 


0.225-0.250 


1.990±0.041 


0.80-0.82 


2.752±0.055 


0.250-0.275 


2.051±0.042 


0.82-0.84 


3.203±0.059 


0.275-0.300 


2.120±0.043 


0.84-0.86 


3.946±0.065 


0.300-0.325 


2.116±0.043 


0.86-0.88 


4.714±0.071 


0.325-0.350 


2.151±0.043 


0.88-0.90 


5.488±0.077 


0.350-0.375 


2.185±0.043 


0.90-0.92 


6.225±0.082 


0.375-0.400 


2.268±0.044 


0.92-0.94 


6.189±0.082 


0.400-0.425 


2.251±0.044 


0.94-0.96 


5.452±0.077 


0.425-0.450 


2.326±0.045 


0.96-0.98 


3.537±0.062 


0.450-0.475 


2.468±0.046 


0.98-1.00 


0.770±0.029 


0.475-0.500 


2.550±0.047 






0.500-0.525 


2.658±0.048 






0.525-0.550 


2.314±0.045 






0.550-0.575 


1.848±0.040 






0.575-0.600 


1.378±0.035 






0.600-0.625 


0.957±0.029 






0.625-0.650 


0.508±0.021 






0.650-0.675 


0.090±0.009 



TABLE I. The measured X3 and X5 distributions with their statistical errors for the three-jet events. The estimated 
systematic uncertainty is 6%. 



cos ^3 


1/NdN/dcoset 


r 


l/NdN/d'<p* 


0.00-0.05 


0.156±0.008 


0.0- 10.0 


0.00726±0.00013 


0.05-0.10 


0.171±0.009 


10.0- 20.0 


0.00780±0.00013 


0.10-0.15 


0.177±0.009 


20.0- 30.0 


0.00780±0.00013 


0.15-0.20 


0.181±0.009 


30.0- 40.0 


0.00673±0.00012 


0.20-0.25 


0.199±0.009 


40.0- 50.0 


0.00561±0.00011 


0.25-0.30 


0.196±0.009 


50.0- 60.0 


0.00440±0.00010 


0.30-0.35 


0.252±0.010 


60.0- 70.0 


0.00376±0.00009 


0.35-0.40 


0.262±0.011 


70.0- 80.0 


0.00327±0.00008 


0.40-0.45 


0.291±0.011 


80.0- 90.0 


0.00325±0.00008 


0.45-0.50 


0.359±0.012 


90.0-100.0 


0.00317±0.00008 


0.50-0.55 


0.416±0.013 


100.0-110.0 


0.00336±0.00009 


0.55-0.60 


0.537±0.015 


110.0-120.0 


0.00384±0.00009 


0.60-0.65 


0.677±0.017 


120.0-130.0 


0.00444±0.00010 


0.65-0.70 


0.928±0.020 


130.0-140.0 


0.00553±0.00011 


0.70-0.75 


1.212±0.023 


140.0-150.0 


0.00686±0.00012 


0.75-0.80 


1.692±0.027 


150.0-160.0 


0.00792±0.00013 


0.80-0.85 


2.160±0.031 


160.0-170.0 


0.00776±0.00013 


0.85-0.90 


2.767±0.035 


170.0-180.0 


0.00722±0.00012 


0.90-0.95 


3.636±0.040 






0.95-1.00 


3.729±0.040 







TABLE IL The measured cos Oz and V'* distributions with their statistical errors for the three-jet events. The estimated 
systematic uncertainty is 6%. 
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^34 


l/NdN/dfis4 


/i35 


l/NdN/dfi35 




l/NdN/dii45 


Q 58—0 60 


0.284±0.018 


0.22—0.24 


099±0 010 


08—0 10 


090±0 010 


0.60—0.62 


0.815±0.030 


0.24—0.26 


0.204±0.015 


0.10-0.12 


0.265±0.017 


0.62—0.64 


1.360±0.038 


0.26—0.28 


0.336±0.019 


0.12—0.14 


0.493±0.023 


0.64—0.66 


2.055±0.047 


0.28—0.30 


0.571±0.025 


0.14—0.16 


0.843±0.030 


0.66—0.68 


2.626±0.053 


0.30—0.32 


0.843±0.030 


0.16—0.18 


1.296±0.037 


0.68—0.70 


3.535±0.062 


0.32—0.34 


1.271±0.037 


0.18—0.20 


1.683±0.043 


70-0 72 


3 719±0 063 


34-0 36 


1 647±0 042 


0.20—0.22 


2.215±0.049 


0.72—0.74 


3.629±0.063 


0.36—0.38 


1.905±0.045 


0.22—0.24 


2.666±0.054 


0.74—0.76 


3.494±0.061 


0.38—0.40 


2.170±0.048 


0.24—0.26 


3.140±0.058 


76-0 78 


3 497±0 061 


0.40—0.42 


2 340±0 050 


26-0 28 


3 317±0 060 


0.78—0.80 


3.491±0.061 


0.42—0.44 


2.465±0.052 


0.28—0.30 


3.684±0.063 


0.80—0.82 


3.514±0.062 


0.44—0.46 


2.628±0.053 


0.30—0.32 


3.716±0.063 


82-0 84 


3 518±0 062 


46-0 48 


2 777±0 055 


32-0 34 


3 592±0 062 


0.84—0.86 


3.625±0.063 


0.48—0.50 


2.914±0.056 


0.34—0.36 


3.434±0.061 


0.86—0.88 


3.467±0.061 


0.50—0.52 


3.297±0.060 


0.36—0.38 


3.293±0.060 


88-0 90 


3 273±0 059 


52-0 54 


3 332±0 060 


38-0 40 


2 916±0 056 


0.90—0.92 


2.296±0.050 


0.54—0.56 


3.750±0.064 


0.40—0.42 


2.582±0.053 


0.92-0.94 


1.267±0.037 


0.56-0.58 


4.024±0.066 


0.42-0.44 


2.369±0.051 


0.94-0.96 


0.482±0.023 


0.58-0.60 


3.919±0.065 


0.44-0.46 


2.102±0.048 






0.60-0.62 


3.390±0.061 


0.46-0.48 


1.887±0.045 






0.62-0.64 


2.679±0.054 


0.48-0.50 


1.584±0.041 






0.64-0.66 


1.945±0.046 


0.50-0.52 


1.176±0.036 






0.66-0.68 


1.136±0.035 


0.52-0.54 


0.907±0.031 






0.68-0.70 


0.301±0.018 


0.54-0.56 
0.56-0.58 


0.565±0.025 
0.154±0.013 



TABLE III. The measured //34, /iss and /i4B distributions with their statistical errors for the three-jet events. The estimated 
systematic uncertainty is 6%. 





l/NdN/dx-i 




l/NdN/dxi 




X 


5 




l/NdN/dx5 


xe> 


l/NdN/dxf. 


0.525- 


-0 


550 


0.14±0.03 





39- 





42 


0.13±0.02 





100- 


-0 


125 


0.15±0.03 





050- 


-0 


075 


0.16±0.03 


0.550- 


-0 


575 


0.43±0.05 





42- 





45 


0.36±0.04 





125- 


-0 


150 


0.30±0.04 





075- 


-0 


100 


0.97±0.07 


0.575- 


-0 


600 


1.12±0.07 





45- 





48 


1.12±0.07 





150- 


-0 


175 


0.46±0.05 





100- 


-0 


125 


1.90±0.10 


0.600- 


-0 


625 


1.51±0.09 





48- 





51 


2.35±0.10 





175- 


-0 


200 


0.81±0.06 





125- 


-0 


150 


3.08±0.12 


0.625- 


-0 


650 


2.18±0.10 





51- 





54 


3.65±0.12 





200- 


-0 


225 


1.10±0.07 





150- 





175 


3.82±0.14 


0.650- 


-0 


675 


2.77±0.12 





54- 





57 


4.49±0.14 





225- 


-0 


250 


1.33±0.08 





175- 


-0 


200 


4.49±0.15 


0.675- 


-0 


700 


2.91±0.12 





57- 





60 


4.57±0.14 





250- 


-0 


275 


1.77±0.09 





200- 


-0 


225 


4.48±0.15 


0.700- 


-0 


725 


3.39±0.13 





60- 





63 


3.94±0.13 





275- 


-0 


300 


2.03±0.10 





225- 


-0 


250 


4.25±0.15 


0.725- 


-0 


750 


3.40±0.13 





63- 





66 


3.25±0.12 





300- 


-0 


325 


2.51±0.11 





250- 





275 


3.69±0.14 


0.750- 


-0 


775 


3.53±0.13 





66- 





69 


2.67±0.10 





325- 


-0 


350 


2.82±0.12 





275- 


-0 


300 


3.26±0.13 


0.775- 


-0 


800 


3.65±0.13 





69- 





72 


2.09±0.09 





350- 


-0 


375 


3.10±0.12 





300- 


-0 


325 


2.95±0.12 


0.800- 


-0 


825 


3.52±0.13 





72- 





75 


1.68±0.08 





375- 


-0 


400 


3.69±0.14 





325- 


-0 


350 


2.25±0.11 


0.825- 


-0 


850 


3.31±0.13 





75- 





78 


1.24±0.07 





400- 


-0 


425 


4.22±0.14 





350- 


-0 


375 


1.82±0.10 


0.850- 


-0 


875 


2.83±0.12 





78- 





81 


0.89±0.06 





425- 


-0 


450 


4.33±0.15 





375- 


-0 


400 


1.35±0.08 


0.875- 


-0 


900 


2.35±0.11 





81- 





84 


0.54±0.05 





450- 


-0 


475 


4.08±0.14 





400- 


-0 


425 


0.83±0.06 


0.900- 


-0 


925 


1.59±0.09 





84- 





87 


0.23±0.03 





475- 


-0 


500 


3.39±0.13 





425- 


-0 


450 


0.52±0.05 


0.925- 


-0 


950 


1.03±0.07 





87- 





90 


0.12±0.02 





500- 


-0 


525 


2.22±0.10 





450- 


-0 


475 


0.15±0.03 


0.950- 


-0 


975 


0.30±0.04 















525- 


-0 


550 


1.08±0.07 

































550- 


-0 


575 


0.52±0.05 

































575- 


-0 


600 


0.11±0.02 













TABLE IV. The measured jet energy fraction distribitions for the four-jet events in their center-of-mass system. Errors are 
statistical only. The estimated systematic uncertainty is 6%. 
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COS 6* 


1 /NdN/d cos ei 


1 /NdN/d cos ex 


1/ NdN/d cos ei 


1/ NdN/d cose; 


0.00-0.05 


0.18±0.02 


0.27±0.03 


0.60±0.04 


0.79±0.04 


0.05-0.10 


0.17±0.02 


0.24±0.02 


0.56±0.04 


0.90±0.05 


0.10-0.15 


0.13±0.02 


0.29±0.03 


0.56±0.04 


1.05±0.05 


0.15-0.20 


0.23±0.02 


0.27±0.03 


0.60±0.04 


0.93±0.05 


0.20-0.25 


0.18±0.02 


0.27±0.03 


0.61±0.04 


0.95±0.05 


0.25-0.30 


0.23±0.02 


0.34±0.03 


0.58±0.04 


0.94±0.05 


0.30-0.35 


0.22±0.02 


0.33±0.03 


0.74±0.04 


0.91±0.05 


0.35-0.40 


0.28±0.03 


0.30±0.03 


0.67±0.04 


0.93±0.05 


0.40-0.45 


0.35±0.03 


0.43±0.03 


0.74±0.04 


0.95±0.05 


0.45-0.50 


0.35±0.03 


0.46±0.03 


0.72±0.04 


0.94±0.05 


0.50-0.55 


0.38±0.03 


0.47±0.03 


0.85±0.05 


1.06±0.05 


0.55-0.60 


0.49±0.03 


0.58±0.04 


0.74±0.04 


1.16±0.05 


0.60-0.65 


0.63±0.04 


0.64±0.04 


0.94±0.05 


1.13±0.05 


0.65-0.70 


0.79±0.04 


0.78±0.04 


0.94±0.05 


1.18±0.05 


0.70-0.75 


0.93±0.05 


0.88±0.05 


1.11±0.05 


1.17±0.05 


0.75-0.80 


1.21±0.05 


1.06±0.05 


1.30±0.06 


1.19±0.05 


0.80-0.85 


1.68±0.06 


1.47±0.06 


1.62±0.06 


1.28±0.06 


0.85-0.90 


2.32±0.08 


2.09±0.07 


2.19±0.07 


1.24±0.06 


0.90-0.95 


3.62±0.09 


3.53±0.09 


2.53±0.08 


0.98±0.05 


0.95-1.00 


5.63±0.12 


5.29±0.11 


1.39±0.06 


0.32±0.03 



TABLE V. The measured jet cosine distributions for the four-jet events in their center-of-mass system. Errors are statistical 
only. The estimated systematic uncertainty is 6%. 



cos UJij 


I /NdN/d cos ujii 


1/ NdN/d cos ^35 


l/NdN/dcosuj36 


l/NdN/dcosco45 


1/ NdN/d cos 0^46 


l/NdN/dcoswid 


-1.0- 


-0.9 


4.557±0.075 


1.180±0.038 


0.550±0.026 


0.043±0.007 


0.407±0.022 


0.782±0.031 


-0.9- 


-0.8 


2.245±0.053 


1.475±0.043 


0.483±0.024 


0.151±0.014 


0.505±0.025 


0.539±0.026 


-0.8- 


-0.7 


1.211±0.039 


1.409±0.042 


0.515±0.025 


0.256±0.018 


0.531±0.026 


0.624±0.028 


-0.7- 


-0.6 


0.815±0.032 


1.150±0.038 


0.577±0.027 


0.300±0.019 


0.535±0.026 


0.596±0.027 


-0.6- 


-0.5 


0.451±0.024 


0.958±0.034 


0.651±0.028 


0.334±0.020 


0.542±0.026 


0.545±0.026 


-0.5- 


-0.4 


0.290±0.019 


0.833±0.032 


0.612±0.028 


0.411±0.023 


0.546±0.026 


0.594±0.027 


-0.4- 


-0.3 


0.197±0.016 


0.774±0.031 


0.651±0.028 


0.422±0.023 


0.566±0.026 


0.603±0.027 


-0.3- 


-0.2 


0.097±0.011 


0.555±0.026 


0.664±0.029 


0.511±0.025 


0.531±0.026 


0.629±0.028 


-0.2- 


-0.1 


0.058±0.008 


0.448±0.024 


0.688±0.029 


0.573±0.027 


0.550±0.026 


0.662±0.029 


-0.1- 


0.0 


0.026±0.006 


0.391±0.022 


0.643±0.028 


0.619±0.028 


0.620±0.028 


0.678±0.029 


0.0- 


0.1 


0.021±0.005 


0.264±0.018 


0.701±0.029 


0.719±0.030 


0.600±0.027 


0.737±0.030 


0.1- 


0.2 


0.022±0.005 


0.208±0.016 


0.697±0.029 


0.745±0.030 


0.636±0.028 


0.738±0.030 


0.2- 


0.3 


0.007±0.003 


0.129±0.013 


0.588±0.027 


0.774±0.031 


0.603±0.027 


0.701±0.029 


0.3- 


0.4 


O.OOliO.OOl 


O.lOOiO.Oll 


0.511±0.025 


0.812±0.032 


0.597±0.027 


0.459±0.024 


0.4- 


0.5 




0.062±0.009 


0.404±0.022 


0.747±0.030 


0.573±0.027 


0.396±0.022 


0.5- 


0.6 




0.037±0.007 


0.360±0.021 


0.745±0.030 


0.513±0.025 


0.311±0.020 


0.6- 


0.7 




0.020±0.005 


0.307±0.019 


0.672±0.029 


0.426±0.023 


0.222±0.017 


0.7- 


0.8 




0.005±0.002 


0.224±0.017 


0.626±0.028 


0.397±0.022 


0.119±0.012 


0.8- 


0.9 




0.002±0.002 


0.144±0.013 


0.420±0.023 


0.262±0.018 


0.054±0.008 


0.9- 


1.0 






0.031±0.006 


0.118±0.012 


0.058±0.008 


0.011±0.004 



TABLE VI. The measured distribution of the cosine of space angles between pairs of jets for the four-jet events in their 
center-of-mass system. Errors are statistical only. The estimated systematic uncertainty is 6%. 
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l/NdN/dlJ,34 




l/NdN/dme 


l/NdN/dfi45 


l/NdN/dfi4G 


l/NdN/dfir^e 


0.00-0.03 














0.03-0.06 








0.02±0.01 


O.OliO.Ol 


0.02±0.01 


0.06-0.09 






0.11±0.02 


0.24±0.03 


0.31±0.04 


0.55±0.05 


0.09-0.12 






0.45±0.04 


0.77±0.06 


0.99±0.06 


1.76±0.09 


0.12-0.15 




O.OliO.Ol 


1.11±0.07 


1.49±0.08 


1.95±0.09 


3.30±0.12 


0.15-0.18 




0.05±0.01 


1.58±0.08 


2.13±0.09 


2.79±0.11 


4.52±0.14 


0.18-0.21 




0.09±0.02 


2.06±0.09 


2.51±0.10 


3.44±0.12 


5.25±0.15 


0.21-0.24 




0.19±0.03 


3.03±0.11 


2.77±0.11 


4.28±0.13 


4.80±0.14 


0.24-0.27 




0.33±0.04 


3.38±0.12 


3.00±0.11 


4.33±0.13 


3.99±0.13 


0.27-0.30 




0.62±0.05 


3.62±0.12 


3.47±0.12 


3.99±0.13 


3.14±0.11 


0.30-0.33 


O.OliO.Ol 


1.23±0.07 


3.90±0.13 


3.36±0.12 


3.61±0.12 


2.32±0.10 


0.33-0.36 


0.06±0.02 


1.57±0.08 


3.49±0.12 


2.94±0.11 


2.88±0.11 


1.65±0.08 


0.36-0.39 


0.09±0.02 


2.14±0.09 


3.07±0.11 


2.99±0.11 


2.13±0.09 


1.05±0.07 


0.39-0.42 


0.12±0.02 


2.91±0.11 


2.82±0.11 


2.72±0.11 


1.40±0.08 


0.60±0.05 


0.42-0.45 


0.25±0.03 


3.53±0.12 


2.03±0.09 


2.26±0.10 


0.84±0.06 


0.29±0.03 


0.45-0.48 


0.61±0.05 


4.35±0.13 


1.59±0.08 


1.47±0.08 


0.37±0.04 


0.09±0.02 


0.48-0.51 


1.06±0.07 


4.92±0.14 


0.81±0.06 


0.95±0.06 


0.02±0.01 




0.51-0.54 


1.64±0.08 


4.60±0.14 


0.25±0.03 


0.26±0.03 






0.54-0.57 


2.74±0.11 


3.63±0.12 


0.04±0.01 


O.OliO.Ol 






0.57-0.60 


3.86±0.13 


2.03±0.09 










0.60-0.63 


3.80±0.13 


0.88±0.06 










0.63-0.66 


4.15±0.13 


0.24±0.03 










0.66-0.69 


3.68±0.12 


O.OliO.Ol 










0.69-0.72 


3.15±0.11 












0.72-0.75 


2.65±0.10 












0.75-0.78 


2.05±0.09 












0.78-0.81 


1.57±0.08 












0.81-0.84 


0.98±0.06 












0.84-0.87 


0.56±0.05 












0.87-0.90 


0.30±0.04 













TABLE VII. The measured distribution of scaled jet pair masses for the four-jet events in their center-of-mass system. 
Errors are statistical only. The estimated systematic uncertainty is 6%. 



Xbz 


1/NdN/dxBz 


cos 


l/NdN/dcos6NR 


0.0-10.0 


0.0073±0.0003 


0.0-0.1 


1.020±0.036 


10.0-20.0 


0.0079±0.0003 


0.1-0.2 


1.068±0.036 


20.0-30.0 


0.0089±0.0003 


0.2-0.3 


0.942±0.034 


30.0-40.0 


0.0102±0.0004 


0.3-0.4 


0.982±0.035 


40.0-50.0 


0.0107±0.0004 


0.4-0.5 


0.993±0.035 


50.0-60.0 


0.0124±0.0004 


0.5-0.6 


0.988±0.035 


60.0-70.0 


0.0132±0.0004 


0.6-0.7 


0.956±0.034 


70.0-80.0 


0.0146±0.0004 


0.7-0.8 


0.986±0.035 


80.0-90.0 


0.0148±0.0004 


0.8-0.9 


0.993±0.035 






0.9-1.0 


1.072±0.036 



TABLE VIII. The measured xbz and cos 9nr distributions with their statistical errors for the four-jet events. The estimated 
systematic uncertainty is 6%. 



25 



